Phase diagrams in nonlocal PNJL models constrained by Lattice 

QCD results 



G. A. Contrera 1 - 2 ' 5 , A. G. Grunfeld 2 ' 3 - 4 , D. B. Blaschke 5 ' 6 
1 Gravitation, Astrophysics and Cosmology Group, FCAyG, UNLP, La Plata, Argentina 
2 CONICET, Rivadavia 1917, 1033 Buenos Aires, Argentina 
3 Department of Physics, Sultan Qaboos University, 
P.O.Box: 36 Al-Khode 123 Muscat, Sultanate of Oman 
4 Departamento de Fisica, Comision Nacional de 
Energia Atomica, (1429) Buenos Aires, Argentina 
5 Institute for Theoretical Physics, University of Wroclaw, 50-204 Wroclaw, Poland 
6 Bogoliubov Laboratory for Theoretical Physics, JINR Dubna, 141980 Dubna, Russia 

(Dated: August 3, 2012) 

We discuss a scheme for extracting the location of the critical endpoint in the 
QCD phase diagram based on nonlocal PNJL models including vector interactions. 
The formfactors of the covariant model are constrained by lattice QCD data for 
the quark propagator. The strength of the vector coupling is adjusted such as to 
reproduce the slope of the pseudocritical temperature of the chiral phase transition 
at low chemical potential extracted recently from lattice QCD simulations. Our 
study supports the existence of a critical endpoint in the QCD phase diagram. For 
the more favorable nonlocal model, with wave function renormalization (set B and 
set C), the critical endpoint is located as (TcePj/^cep) = (129.8 MeV, 276.6 MeV) 
and (Tcep> A'cep) = (69.9 MeV, 319.1 MeV) for set B and set C, respectively. 

The conjecture for the existence of a critical endpoint (CEB) of first order phase tran- 
sitions in the QCD phase diagram is the basis for recent as well as future beam energy 
scan (BES) programs in relativistic heavy-ion collision experiments at RHIC, SPS, NICA 
and FAIR which try to identify the parameters of its position (Tcep, A*cep)- The theoretical 
situation is very unsatisfactory since the predictions for this position form merely a skymap 
in the T — \x plane [1] . 

The now well-established results from lattice QCD that at zero and small chemical poten- 
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tial /i, predict coincident chiral and deconfinement crossover transitions at a pseudocritical 
temperature of T c (0) = 154 ± 9 MeV for 2+1 flavours [2] and a value of T c (0) ~ 170 MeV 
for two flavours [3|. Then, the underlying idea is to use effective theories based on QCD, 
which reproduce lattice results at vanishing /i, and see whether the predictions can be ex- 
trapolated at high chemical potential. That leaves us with a variety of possibilities for the 
phase structure at non zero fi, among them we want to mention the following ones 

• no CEP at all since the transition is crossover in the whole phase diagram, 

• no CEP, but a Lifshitz point [3], 

• one CEP, but with largely differing predictions of its position PQ, 

• second CEP [6HH], 

• CEP and triple point, possibly coincident, due to another phase (i.e. colour supercon- 
ducting [H] or quarkyonic [TO] matter) at low temperatures and high densities. 

This situation is rather unsatisfactory in view of the upcoming experimental programmes. 
In the present work we discuss the existence and location of a CEP within the frame of 
nonlocal chiral quark models coupled to the Polyakov loop potential on the selfconsistent 
meanfield level. 

The Lagrangian for the SU(2) f non local models, including vector channel interactions, 
is given by 

£ = q(ilp-m )q + £ int + U($) , (1) 

where q is the Nf = 2 fermion doublet q = (u,d) T , and m is the current quark mass (we 
consider isospin symmetry, that is mo = m u = m^). The covariant derivative is defined as 
Dp = dp — iAp, where Ap are colour gauge fields. 
The nonlocal interaction channels are given by 



int ~ 2 



i G 

ja(x)ja(x) ~ jp(x)jp(x) ^jv(x)j V (x), (2) 
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where the nonlocal currents are 

z 



Ja{x) = J d 4 z g(z) q(x+^j T a q \x - - 

3p( x ) = / d4z f( z ) Q ( x + |) i ( 3 

J p 

J d 4 z g(z) ?(x + ^) 7 %(x- |) (3) 



z 
2 



3v{x) 

with T a = (r5,Tp) = (1,z7 5 t) for scalar and pseudoscalar currents respectively, and 

^ y j — 

u(x') d v(x) = u(x')d x v(x) — d x iu(x')v(x). The functions g(z) and f(z) in (K3) are nonlocal 
covariant form factors characterizing the corresponding interactions. The scalar-isoscalar 
component of the j a {x) current will generate the momentum dependent quark mass in the 
quark propagator, while the "momentum" current, jp(x), will be responsible for a momen- 
tum dependent wave function renormalization (WFR) of this propagator. Note that the 
relative strength between both interaction terms is controlled by the mass parameter k v 
introduced in Finally, jv(x) represents the vector channel interaction current, whose 
coupling constant Gy is usually taken as a free parameter. 

In what follows it is convenient to Fourier transform into momentum space. Since we are 
interested in studying the characteristics of the chiral phase transition we have to extend the 
effective action to finite temperature T and quark chemical potential fi [33J. In the present 
work this is performed by using the Matsubara and imaginary time formalisms. Concerning 
the gluon fields we use the same prescription as in previous works [TU H2] , but in our present 
case we have chosen a //-dependent logarithmic effective potential described in [13J. 

T, n) = (a T 4 + a ly u 4 + a 2 TV )$ 2 + a 3 T 4 In (1 - 6$ 2 + 8$ 3 - 3$ 4 ), (4) 

where ao = —1.85, ai = — 1.44xl0~ 3 , a 2 = —0.08, a 3 = —0.40. In the present work we set 
the To parameter by using the value corresponding to two flavours To = 208 MeV, as it has 
been suggested in [TJ], and used in subsequent approachs, including the nonlocal PNJL 
[To] and Polyakov loop-DSE models [16]. 

Finally, to fully specify the nonlocal models under consideration we fix the model pa- 
rameters as well as the form factors g(q) and f(q) following [121 E]> i- e - considering two 
different types of functional dependencies for these form factors: exponential forms (set A 
and set B) and Lorentzians (Set C). In order to compare with some recent results |4J, we 
have included a local form factor, using for this case the parameters in [18J. 
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The model inputs can be constrained with recent results from lattice QCD studies. In 
particular, the form factors of the nonlocal interaction [17J have been chosen such as to 
reproduce the dynamical mass function M(p) and the WFR Z(p) of the quark propagator 
in the vacuum |19j . In Figure [l] we show the shapes of normalized dynamical masses and 
WFR for the models under discussion here, i.e., the nonlocal models of set A (rank-one), set 
B and set C (rank- two) as well as the local limit. This figure generalizes the corresponding 
Figure [l] in [17] . by including the local limit case and showing the line Z(p) = 1 for both, 
set A and the local case. It is easy to recognize the better agreement between the lattice 
results and the more complete model, namely nonlocal PNJL with WFR, in its two form 
factor parametrizations given by set B and set C. 
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FIG. 1: Normalized dynamical masses and wave function renormalization for the different form 
factors under study, fitted to lattice data [19], extended from [T7] by inclusion of local model [18J. 
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Within this framework the mean field thermodynamical potential f2 MFA results 
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MFA 
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4T EE 



°l + «p ^2 

2Gc 



d 3 p 
(2tt) 3 



2Gv 



In 



(P^) 2 + M 2 (p^) 



■W($,T). 



(5) 



where M(p) and are given by 

M(p) 
Z(p) 



Z{p) [m + ai g(p)] 
[1-^2 /(p)]" 1 . 



In addition, as in 



we have considered 

2 



(2n + l)wT - i/x + 0, 
= diag( 



where the quantities C are given by the relation 

c 



r ,(f) g ,(j) b ). Namely, 



(6) 
(7) 

c 03 with 



1, —1,0 for r,g,b, respectively. In the case of (p^pj we have used the same definition 
as in ([7]) but shifting the chemical potential according to 

fi = n -wg(p) Z{p). (8) 

The f2 MFA turns out to be divergent and, thus, needs to be regularized. For this purpose 
we use the same prescription as in [121 120] • The mean field values cxi^, u and 03 at a given 
temperature or chemical potential, are obtained from a set of four coupled "gap" equations 
which come from the minimization of the regularized thermodynamical potential, that is 



(9QMFA 



reg 



uiL reg 



0. 



(9) 



dui da 2 du dcj) 3 

We also want to include in our analysis the results arising from a local PNJL model. 
For that purpose, we started from the Lagrangian in [2T] with two flavours instead of three. 
Moreover, we consider that the chemical potential is shifted by 



H = yU —u. 



(10) 

We use the same set of parameters as in [T8]. 

The vector coupling constant Gy is considered a free parameter which in the mean field 
approximation (MFA) may be adjusted such as to reproduce the behaviour of TJjjl) which 
has recently been obtained by Taylor expansion technique in lattice QCD [22] 



W/T c (0) 



K^/T c ) 2 + 0[iti/T c y 



(11) 
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with k = 0.059(2) (4) being the curvature. 

Note that this result is not yet based on continuum extrapolated lattice results, so that 
discretization errors has to be expected. However, we are interested to present here a 
scheme for constraining effective QCD models. The quality of predictions can subsequently 
be improved by using better lattice QCD data and constant discretization schemes. 

To determine the curvatures we plotted the pseudocritical temperatures of the crossover 
transitions at low values of \xjT for different ratios 7]y = Gy /Gg. An example of this is shown 
in Figure [2] for set C (the corresponding plots for the other models are qualitatively very 



similar). The fit of the lattice QCD results (11) is also shown. The grey zone corresponds 



to the error in the coefficient k obtained in [221. 
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FIG. 2: Chiral crossover transitions at low values of fj,/T for different values of strengths of the 
vector coupling rjy = Gy/Gs for Set C. The dashed line corresponds to the lattice QCD prediction 
of k = 0.059(2)(4) [22]. 



In Figure [3] we compare the lattice result with the values for the coefficient k obtained 
within the nonlocal PNJL models and the local one. There, the horizontal line corresponds 
to the lattice QCD prediction of k = 0.059(2) (4) [22J and the grey zone represents its error. 
Note that for the more complete model (set B and set C) the curvatures are closer to each 
other than in the case of set A and local ones. 

We would like to remark that while the analysis of k has been performed for Nf = 2 + 1 
simulations, the chemical potential /x concerns only the two light flavours. Therefore, our 
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extraction of k from the nonlocal PNJL models for the 2-flavour case may be in order. 
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FIG. 3: Curvature k of the pseudocritical temperature T c (fi) of the chiral crossover transition at 
low values of \xjT . The horizontal line corresponds to the lattice QCD prediction of n = 0.059(2) (4) 

In Figure [4] we compare the two different classes of nonlocal models discussed in this work. 
The nonlocal model for set B and set C cases, contains WFR and dynamical (running) quark 
mass effects and thus is closer to full QCD. Therefore, we suggest that statements about 
the existence and location of the CEP within set A and the local model should be less 
trustworthy than those of set B and set C. 

In Figure [4] we have also shown a grey zone corresponding to the range of r\y values 
obtained by taking into account the error in the lattice determination of k [22J. The error 
bars in the CEP's indicate the distances to the CEP positions for the corresponding r\y 
values that fit the error limits. Note that in the local model, for fitting the lattice QCD 
value a larger vector coupling is required than in the nonlocal ones. The consequences for 
predictions about the existence and location of the CEP are dramatic, see Figure [5} On 
the other hand, we found that using the Polyakov loop potential Q from [13J, a crossover 
region and a CEP can be obtained for set A, even for To = 208 MeV, contrarily to what has 
been reported in [T5l H~6] . where the Polyakov loop potential from [23] has been used |34j. 

The results summarized in Figure [5] indicate that the absolute value of the critical tern- 
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FIG. 4: Phase diagrams for the non local models for the values of rjv that best fit the lattice QCD 
results in [22]. 

perature T c (0) of nonlocal covariant PNJL models is rather insensitive to the choice of the 
form factors parametrizing the momentum dependence of dynamical (running) mass func- 
tion and WFR of the quark propagator as measured on the lattice at zero temperature pjj] , 
whereas the position of the CEP and critical chemical potential at T = strongly depends 
on it. Nevertheless, the value of T c (0) in the local model is significantly different (larger) 
than in the nonlocal ones. 

The region of CEP's suggested by our study is between the results for set B and set C, 
i.e. (Tcep^cep) = (129.8 MeV, 276.6 MeV) and (69.9 MeV, 319.1 MeV), respectively. 

This means that the search for CEP signatures in the BES programs is justified and 
should be continued. The energy range of the NICA and FAIR facilities shall be particularly 
promising. As a characteristic feature of the region around the CEP any phase transition 
signals based on changes in the bulk properties are not pronounced. Therefore, sensible 
observables are required for the identification of the CEP position like, e.g., higher moments 
of the identified particle distributions [MJ |25] or the flow of nuclear clusters [26] . 

In conclusion, we have presented progress in the extraction of the CEP location in the 
QCD phase diagram based on nonlocal PNJL models constrained by lattice QCD. Note 
that in this exploratory study the lattice QCD data for the quark propagator data and 
for the curvature of the pseudocritical line are obtained with different lattice actions. A 
more consistent study should be based on the same discretization of the action, provide a 
continuum extrapolation and work with physical quark masses. 
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FIG. 5: (Colour online) Phase diagrams with (pseudo)critical temperatures T c (p) and critical points 
for nonlocal PNJL models (set A-C) compared to the local one. Dashed (full) lines correspond to 
crossover (first order) transitions. The vector coupling strength r\y is chosen so that these models 
reproduce the lattice QCD result k = 0.059(2) (4) |22j for the curvature at low ^-values. The 
corresponding values for T c (/x = 0) (in MeV) are 169.9, 171.3, 173.2 and 200.9, for Sets A, B, C 
and local, respectively. The highlighted region denotes the CEP position favored by the present 
study. 

As a next step it is necessary to investigate the robustness of the results of the nonlocal 
PNJL models when modifying the choice of the Polyakov-loop potential taking into account 
recent developments [2TH29"] and, in particular, when going beyond the mean field approx- 
imation. A scheme for going beyond the meanfield in nonlocal PNJL models by including 
hadronic correlations (bound states and their dissociation in the continuum of scattering 
states) has recently been developed [30] and shall be generalized for studies of the chiral and 
deconfinement phase transition in the QCD phase diagram. A key quantity for such studies 
will be the hadronic spectral function and first results using a generic ansatz [3T] for joining 
the hadron resonance gas and PNJL approaches are promising [32 J. 
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